Noise-Induced Phase Separation: Mean— Field Results 
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I. INTRODUCTION 

Many theoretical and experimental observations con- 
firm nowadays the existence of noise-induced order. Phe- 
nomena such as noise-induced transitions [Q, stochas- 
tic resonance and noise-induced transport ||^ are ex- 
amples of the ordering features of fluctuations in purely 
temporal dynamical systems. Additionally, recent years 
have witnessed an increasing interest on noise-induced 
phenomena in spatially-extended systems (see Q for a 
recent review). Some of the topics studied in this re- 
spect include noise-induced patterns [^||^] , noise-induced 
phase transitions spatiotemporal stochastic reso- 

nance [|lO|,0 , noise-induced fronts noise-supported 
traveling structures in excitable media p3| and noise sus- 
tained convective structures ||l^,|l5| . We are concerned in 
this paper with the phenomenon of noise-induced phase 
separation, recently observed in systems with conserved 
dynamics JTot . 

Several analytical methods have been used so far Q 
to examine the above-mentioned spatiotemporal prob- 
lems. By way of example, the stability of a homoge- 
neous state with respect to small perturbations of arbi- 
trary wavenumber can be analyzed in a linear approx- 
imation. Such a linear stability analysis shows that 
pattern-forming transitions are nontrivially affected by 
multiplicative noise From a more fundamen- 

tal point of view, systems exhibiting phase transitions 
in a statistical-mechanics sense can be investigated by 
means of the dynamic renormalization group |]l8| , |l9| , 
which shows that under certain conditions a new gen- 
uine nonequilibrium universality class arises due to the 
presence of multiplicative noise [^,^ . A third fruitful 
approach is based on the well-known mean-field approx- 
imation, widely used in the context of equilibrium sta- 
tistical mechanics, and that has been recently extended 



to nonequilibrium systems under the influence of exter- 
nal noise 0,^. In this context, the approximation is 
introduced by assuming that the interaction between a 
certain spatial point and its neighbors occurs through a 
mean value of the field, which corresponds to its statis- 
tical average at the given point. This approach has led 
to the prediction of noise-induced ordering and disor- 
dering phase transitions (NIOTs and NIDTs), which has 
been successfully verified (at least qualitati vely) by nu- 
merical simulations in different models |0, 



21 



The 

advantages of this procedure as compared to, e.g., linear- 
stability approaches, lie on its ability to describe the sys- 
tem arbitrarily far from the transition point and to take 
into account the influence of spatial coupling strength, 
that arises naturally in discretized systems. In this way, 
the mean-field analysis can successfully explain the exis- 
tence of successive NIOTs and NIDTs (also called reen- 
trant transitions in this context) as a single control pa- 
rameter is varied. 

The aim of this paper is to perform a somewhat de- 
tailed study, using the mean-field approximation tech- 
nique, of the phenomenon of noise-induced phase sepa- 
ration. This phenomenon has been recently predicted by 
a linear stability approach and confirmed by numerical 
simulations Qjl^. It arises in spatiotemporal systems 
whose dynamics is conserved, in the sense that the spa- 
tial average of the field does not vary with time, but de- 
pends only on the initial conditions of the system. Due 
to this fact, a standard mean-field approach cannot be 
applied in this case, because no change in the mean field 
will be observed as a given control parameter is varied 
(and hence no phase transition can be found in this way) . 
Therefore, an extension of the procedure is needed in or- 
der to handle this situation. The present work is devoted 
to developing such an extension, and applying the re- 
sults to the particular case of noise-induced phase sep- 
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aration mentioned above. The outline of the rest of the 
paper is the following: Section II introduces the general 
system that will be investigated, along with the particu- 
lar model to which the obtained results will be applied. 
A comparison between conserved and nonconserved dy- 
namics is also briefly sketched. Section III reviews the 
mean- field procedure for nonconserved systems, and ex- 
tends it to include the effect of spatial correlation of 
the external noise. Section IV introduces the general- 
ized mean-field approach for conserved systems. Section 
V discusses the limit of strong spatial coupling of the 
procedure, and compares the corresponding results with 
those coming from linear stability analysis. Throughout 
all these Sections, a comparison with respective numeri- 
cal simulations of the complete model is made. Finally, 
some conclusions are stated in Section VI. 



components may undergo, for some values of the con- 
trol parameters, a separation process which, depending 
on the initial relative concentrations of each component, 
takes the form of spinodal decomposition or nucleation 
psf . In this paper we will be mainly concerned with 
the conserved case, although a comparison will also be 
made with the corresponding nonconserved case, defined 
by r = 1 (known as model A ). 

Even though the theoretical approach that will be pre- 
sented here is quite general, our results will be applied, 
for the sake of clarity, to the particular Ginzburg-Landau 
form of the free energy F ^ 



F ^ dx 



n0) + 2iv0p 



(4) 



where the local potential V{(t>) is 



II. CONSERVED AND NONCONSERVED 
DYNAMICAL MODELS 

The spatiotemporal dynamics of a nonequilibrium sys- 
tem in the presence of both internal and external noise 
can be described by the following Langevin equation ||2^ ] 
for the time evolution of the d-dimensional scalar field 
t): 



d(j){x, t) 
dt 



6F 



+ Vix,t), (1) 



where the additive noise 77(2?, t) is Gaussian and white, 
with zero mean and correlation 



{ri{x, t) v{x',t')) = 2 £ r (5(f - x') 5{t - t') 



(2) 



The intensity of the noise is measured by the parame- 
ter £. The existence of the factor F in correlation is 
a sign of the internal character of this noise, in whose 
only presence (g — 0) the system can exhibit equilibrium 
properties. The multiplicative noise term ^(x, t), on the 
other hand, is external and brings the system out of equi- 
librium. It may arise, for instance, from a fluctuating 
control parameter. It is also Gaussian with zero mean, 
but its correlation will be assumed in principle to have a 
nontrivial structure in space: 



(^(f , t) ^x', t')) = 2 a' c(|f ~ x'\) 6{t - t') , 



(3) 



1/(0) 



(5) 



In the absence of noise sources, the behavior of this po- 
tential is the following: for a < the homogeneous trivial 
solution (f) — Q \s stable, whereas for a > that solution 
becomes unstable. This instability gives rise either to 
a phase transition towards an ordered (ferromagnetic) 
phase in the nonconserved case, or to a phase separation 
process in the conserved case. 

The external noise will be taken to be coupled to the 
field according to 



(6) 



which corresponds to allow the control parameter a in (|^) 
to fluctuate in space and time. We will use the following 
Gaussian spatial correlation function 



exp 



2A2 



(7) 



whose width A characterizes the correlation length of the 
noise. The normalization is such that in the limit A — > 
this correlation goes to a delta function and ^(x, t) be- 
comes a spatial white noise with intensity . 



III. MEAN-FIELD APPROACH FOR 
NONCONSERVED DYNAMICS 



where c(|a; — x'|) is the spatial correlation function of the 
external noise and is its intensity. 

Different and physically motivated choices for F will 
lead to a variety of dynamical and steady state phe- 
nomenologies. The particular case of F = — (called 
model B in the literature of critical phenomena) is appro- 
priate to describe a system in which the global quantity 
J 4){x, t)d'^x is conserved in time. Physical realizations 
of this system include the case of phase separation in bi- 
nary alloys. In this case, an initial mixture of the two 



We now review the main points of the mean-field 
approach in its application to nonconserved order- 
parameter systems (model A), in order to clarify the ex- 
tension to conserved dynamics that will be presented in 
the next Section. We begin by discretizing the field equa- 
tion (§) with F = 1 in a regular c?-dimensional lattice of 
mesh size Ax = 1, and lattice points afi, . . . , xm 



D ~ 



^J+^l^{t) +g{^i)Ut), (8) 
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where = 4'i^i)i fi'Pi) = ~V{4'i)j and only one index 
is used to label the cells, independently of the dimension 
of the lattice. Dij accounts for the discretized Laplacian 
operator 



nn(i),j 



(9) 



where nn{i) represents the set of all the sites which are 
nearest neighbors of site i. 

The discrete noises r]i{t) and ^i{t) are still Gaussian 
with zero mean and their correlations are given by 



and 



where 



(10) 



(11) 



is a convenient discretization of the function 



c{\x — x'\) and specific values such as cq, ci have to be 
computed numerically UJl^ when needed. For the white 
noise case (A = 0) and with the mesh size chosen Arc = 1 
it is Co = 1, ci = 0, whereas for large A, cq scales roughly 
as Co oc A~'^. 

The corresponding Fokker-Planck equation, in the 
Stratonovich interpretation, for the multivariate proba- 
bility density Pi^uh, ■ ■ ■ ,t) = P{{<p}, t) is [| 



dP 

'dt 



E 



d 



— e - 



oE 



d 



P. 



(12) 



In order to get the evolution equation for the single- 
site probability distribution P{(f)i,t), defined as 



P(0»,t) 



(13) 



we integrate Eq. (^^ over all the variables except (j^i. 
Vanishing of the probability for the field going to ±00 
leads to 



^ [5(0,)P(W,t)]-O j^t (14) 



and using the standard definition of the conditional prob- 
ability, one gets 



E 

j£nn{i) 



Y]^ d(f)k 



k^i 



J2 J dq^, Piq^,\<P,,t) 



2d (0(0)0. P(0.,i) (15) 



which defines {(f>{t))^- as a nearest-neighbor conditional 
average. Thus we finally find that the one point steady 
probability distribution follows the simpler but still exact 
equation. 



^P{^^,t) 

dt 



-^[/(0,O + ^((0(0)0.-0.) 



■9{4>i) 



P{(l),,t). (16) 



The mean-field approximation consists in assuming 
that the conditional average in the last equation is re- 
placed by M 



(17) 



which is equivalent to doing directly the following as- 
sumption at the level of the Langevin Eq. (H), 



A,0,(i)-(0»(t))-<^,:(i) 



(18) 



Using this approximation, and imposing the condition of 
stationary probability distribution with no flux, we get 
that the single-site steady distribution satisfies, 



/(</>) -|-i?((0)st - 0) - 



d 



^^st(0)=O, 



(19) 



where subscript i has been dropped for simplicity. 

The solution of the previous equation can be easily 
written down: 



Pst(0, (0)st) = iVexp 



1 



0-2 Co g'^{4>') + e 
[/(0') + D ((0>st - cp') - CO g{cl>') g'{^')] , (20) 

where N is an appropriate normalization constant. The 
above solution is only formal, because {(f))st depends on 
the probability distribution itself. However, both Pgt and 
{(j))st can be determined by means of the following self- 
consistency relation, which is a signature of the mean- 
field approach. 



(j)P,t{(l>,{cj))st)d<i). 



(21) 



We now apply our results to the particular model de- 
fined by (^-(0)- The solution of the self-consistency 
equation (2l|) in this case is plotted in Fig. |^ versus the 
control parameter a for three different values of the noise 
correlation length A. Note that the existence of multi- 
plicative noise shifts the critical point towards negative 
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values of a, which indicates the ordering character of the 
external noise. This noise-induced phase transition has 
been substantially studied in the past in the case of white 

Figure also displays the re- 
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external noise 
suits obtained by numerical simulations of the complete 
model (|^) for a two-dimensional square lattice and us- 
ing the mean-field results as the initial conditions. It 
can be seen that mean-field results give the correct qual- 
itative behavior of the system, and are also quantita- 
tively right far from the critical point. The agreement 
between mean-field predictions and simulations close to 
the critical point improves when the correlation length A 
increases. 

Notice that, in this mean-field approximation, the only 
effect of the finite correlation length A of the noise shows 
up in the value cq < 1 (see Eq. (po|)) of the correla- 
tion function at zero distance. As mentioned before, this 
value decreases with increasing A. In other words, for 
nonconserved dynamics the disordering effect of the spa- 
tial correlation of multiplicative noise in the mean-field 
approximation arises only through a decrease of the ef- 
fective noise intensity. 



IV. MEAN-FIELD APPROACH FOR 
CONSERVED DYNAMICS 

The mean-field approach discussed above cannot 
be straightforwardly extended to deal with conserved- 
order-parameters systems, because in these cases the 
mean field {(j)) is constant in time, depending only on 
the initial conditions and not on the control parameter a. 
We now introduce a generalized mean-field approxima- 
tion that overcomes such a restriction. The main ideas 
underlying this extension will be first presented in the 
deterministic model B. 



A. Deterministic dynamics 

In the absence of all noise sources, model B takes the 
form 



dt 56 



(22) 



This model evolves in time under the following restriction 



1 

V 



dx 4>{x, t) = (f>Q , 



(23) 



being fixed by the initial conditions. The phenomenol- 
ogy of this model is well known : there is a transition 
point aT{4>a), such that for a < ariipo) the homogeneous 
state (j) = 4)0 is stable, whereas for a > aT{4>o) the sys- 
tem separates in two bulk phases, ^^nd (f>2, fulfilling 
that the spatial average of (j) is also equal to 0o. The 
transition from an homogeneous state to a two-phase 



state is critical {i.e., of second order) for (/iq = 0, so that 
aT(<^o = 0) = flc. 

In order to determine both axicfio) and ac, we look 
for the steady-state solutions of Eq. (p2|). These solu- 
tions fulfill the Laplace equation V^|^=: 0. The an- 
alytical and bounded solution is ^ = h, being h a 
constant. Therefore, the steady states of model B can 
be interpreted as the minima of an effective potential 
Foff — F — h J dx(j), and coincide with the steady states 
of model A with an external control field h. Following 
reference [ p7[ we call h the constant effective field of the 
system. For equilibrium systems h is nothing but the 
chemical potential. Moreover, h is not an arbitrary con- 
stant, and it has to be determined by imposing the con- 
servation law, Eq. (23). Substitution of the Ginzburg- 
Landau form in the discretized version of h — leads 
to 



-a(f)i + 0; 



D 
2d 



(24) 



We now need to consider separately the subthreshold 
and superthreshold situations: 

• In the subthreshold (homogeneous) case the condi- 
tion = 00: Vi, has to be verified, and therefore 
Eq. (|3) reduces to 



h = -acpo + 



(25) 



Hence the value of h does depend on the initial 
condition in the subthreshold situation. 

• Above the transition point (not yet determined), 
the steady state of the system is not globally homo- 
geneous, since the field separates in two bulk phases 
with values 0i and (j)2, respectively. The fraction 
X of system in phase (jji is given by the lever rule: 
X(j)i -|- (1 — x)(l>2 = </'o- For a general free energy, a 
Maxwell-type construction would give us the value 
of h. In the case of a locally symmetric free-energy 
(such as the one defined in Eqs. (|4|)) and (1), a sim- 
pler argument can be used: each phase has to sat- 
isfy Eq. (^) with (jji equal to the field value of the 
corresponding phase, either 0i or (j)2, and, since by 
the symmetry of the free energy, these two quanti- 
ties verify (f>i = —4>2, h must be zero. Consequently 
we get 



(26) 



which are the solutions of the deterministic model 
A for a value of the external control field h ^ 0. 

Just at the transition point, there is a unique phase 
= 00 and h is identically zero. Thus the transition 
line (also called in this context coexistence line) is given 
by 



(27) 
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We also note that for <j>Q = 0, the critical point is ob- 
tained: Ut = CLc = 0. 

We will now show that the concept of the constant 
effective field can be used to generalize the mean-field 
approximation to conserved systems with noise. 



B. Noise— induced phase separation 

We now add stochastic sources to model (p^), in the 
form of both an internal additive noise and an external 
multiplicative one. The resulting model is represented by 
Eqs. (0)~(|) with r = -V^ The discretized version of 
this model is 



d(t)i 
~dt 



] \ k 



(28) 



with f{(j)j) = —V'{(pj), as before. The correlation of the 
additive noise is now 



(29) 



and that of the multiplicative noise was already intro- 
duced in Eq. (|ll|). The corresponding Fokker-Planck 
equation, in the Stratonovich interpretation, for the mul- 
tivariate probability density P{{(j)},t) is in this case 



dP 
'dt 



i,j ^ \ k 



" ^ 9^ + ^^^'^ ^ W "l^-'-l ^^'^'■^ ^ ■ 

r,s / 

As done in the nonconserved case, we now integrate 
Eq. (^) over all the variables except i^i, in order to get 
the evolution equation of the single-site probability dis- 
tribution Pi(t)i,t) (see Eq. (|l|)) 



dPiq}„t) d 



dt 



— J2D.,,{M,)^^P{^.,t), (31) 



where 



D 



d_ 

d(t)j 



+ '^^9{(l>j) E -Qr^'^^^\]-r\9{(t>r). (32) 

If we impose the condition of stationarity probability dis- 
tribution with no flux, {Mj)^^ must satisfy 



(33) 



In the deterministic conserved case it has been shown 
that the solution of this equation is the constant effective 



field, {Mj}^. = —h. We can now take j = i and perform 
the conditional average of Mi. If we consider the expres- 
sion analogous to Eq. ( [T^ for the multiplicative noise 
term and make the standard mean field approximation 
( [TtI ) we arrive at 



hPstW = 
f2dCT^g((/>) 



ci5((0>st)^-co^5(0) 



d 



Pst{(t>), (34) 



where subscripts have been dropped again for simplic- 
ity. In the derivation of this equation a generalization of 
the mean-field approximation for the nearest-neighbor 
conditional average of function gicj)) has been applied, 
namely: 



(35) 



In principle this is an uncontrolled approximation whose 
validity needs to be assessed by the numerical simulations 
which will be presented in what follows. 

The solution of Eq. (^4|) yields the following stationary 
probability distribution 

P,t(0, {4')st:h) = 



A^exp 



2da^g{cl>'){cog{cl>')-cigmst))+ e 
[/((/.') + D ((0),t - 0') - 2da2co<?((/.')<?'(0') + h] , (36) 

where h and ('/>)st are parameters to be determined sclf- 
consistently. 

We now particularize the result obtained above to the 
Ginzburg-Landau model defined by Eqs. (|^)-(0). In this 
case, the stationary single-site probability distribution is 

P,t(</', {^)st,h)=Nx 



exp 



/ 



D - 2da^co) (j)' - (j)'^ + D{(l3)st + h 
2da2 (co</)'2-ci(0),t 0')+e 



d<j)' . (37) 



We now need to determine the unknown constants h 
and (0)st- Similarly to the deterministic case studied 
above, we consider separately the subthreshold and su- 
perthreshold situations. We recall at this point that the 
mean-field approach presented above is local, and leads 
to an expression for the probability distribution of the 
field at a given site of the lattice as a function of h and 
of the mean field (0)st in the neighborhood of the given 
cell. In the homogeneous case (a < ax), this mean field 
is the same everywhere, and it is equal to (j^o- Hence only 
h is left to be evaluated, what can be done by means of 
the generalized self-consistency relation 
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(38) 



with {(j))st = 4>a- 

For a > ax the system has two phases, and thus there 
are two different local mean values, corresponding to each 
of the bulk phases {{41)1 and (0)2)- Because of the sym- 
metry of -F and g, these values must satisfy = —('/') 2- 
Therefore, since h needs to be the same for the two 
phases, and given the form of the numerator in (^^, h 
must be zero in this ordered state. Hence only the val- 
ues of the local (symmetrical) mean fields and (0)2 
need to be determined for a > ut- This can be done as 
in the case of model A, solving the self-consistency rela- 
tion (|2^) using the steady probability given in (^^ with 
h = 0. For nonsymmetric functions F{(j)) and g{(t>) a pos- 
sible extension of the Maxwell rule is to choose h in such 
a way that the two solutions of Eq.(^8|) have the same 
value for the stationary probability -Pst((^)stj (0)sti h). 

The bifurcation diagram resulting from the application 
of the self-consistency relations is plotted in Fig. ^ for 
three different values of the multiplicative-noise correla- 
tion length. Numerical simulation results of the complete 
model (28) are also shown. Mean field results have been 
used as the initial conditions of the numerical simulations 
letting each phase evolve until its stationary value. The 
effects of the intensity and correlation length of the mul- 
tiplicative noise are qualitatively the same as in model 
A: the noise-induced shift of the transition point, in the 
direction of enhancing order in the system, increases with 
noise intensity and decreases with correlation length. 

Figure || shows the values of the constant effective field 
h obtained numerically by imposing the self-consistency 
relation ( ^8|) until h vanishes, for a nonzero initial con- 
centration (pQ = 0.2. The corresponding value of the 
control parameter a at which h first becomes zero is the 
transition point ax- Results have been plotted for the 
deterministic case (which can be calculated analytically, 
see Eq. (p5|)), the case with just additive noise, and three 
cases with also multiplicative noise for different correla- 
tion lengths, corresponding to the situations shown in 
Fig. ||. The noise-induced shift of the transition point 
and the infiuence of the noise correlation length as well 
as the disordering role of the additive noise (reflected in 
the shift of the transition point towards the right when 
only additive noise is considered - dotted line) can be 
clearly seen. 

A comment on the comparative influence of multiplica- 
tive noise on conserved and nonconserved dynamics is 
worth to be made at this point. We note that, in the 
ordered state {h — 0), the single-site probability distri- 
bution of the conserved model ( ^6| ) in the presence of 
white multiplicative noise (cq = 1 and ci = 0) reduces to 



P,t(0, (0)st) = A^exp 



1 



which should be compared with the corresponding ex- 
pression ( po|) for the nonconserved case with cq = 1. One 
can easily see that multiplicative noise has a stronger ef- 
fect on the conserved model than on the nonconserved 
one, since in the former case the noise intensity is multi- 
plied by a factor 2d. In the particular case in which the 
two noise intensities of the nonconserved (A) and con- 
served {B) models are related by 



2da|, 



(40) 



2dcr2 g2(0') +e 
[/(^') + D ((0),t " 0') - 2da^g{cp')g'{(b')] , (39) 



the two models are equivalent above the transition 
point. However, this equivalence disappears in the 
case of colored multiplicative noise, because of the term 
2da'^ cig{(j)')g({(j))st) appearing in Eq. (^). This differ- 
ent dependence indicates that spatial correlation of the 
noise is more relevant for the conserved model than for 
the nonconserved one, where correlation length of the 
noise produces only a shift of the transition point [ p^ . 
A comparison between the results of model A and B is 
shown in Fig. ^, for both A = and A ^ 0. Noise in- 
tensities have been chosen here to verify expression (|4^), 
so that in the white-noise case, mean-field results coin- 
cide for the two models. Mean-field results are in better 
agreement with simulations in the case of model B. 

Finally, we now address the issue of whether a reen- 
trant noise-induced phase transition towards disorder 
arises in the conserved model B. Previous works have 
shown the existence of such a transition for nonconserved 
models This means that for fixed values of a and 

D, when increasing the multiplicative noise intensity, the 
system goes first through a phase transition from disor- 
der to order (NIOT) and then, for higher values of the 
noise, it experiments another transition back to disorder 
(NIDT). These two transitions can only be found when 
increasing , instead of a or Z?. 

Mean-field theory predicts the existence of reentrant 
transitions also for model B, as shown in Fig. |[ This 
figure shows the behavior of the mean field {4>)st versus 
multiplicative noise intensity with two different correla- 
tion lengths (A = and A = 0.5) for model A and B. 
Again, the noise intensities for the two models have been 
chosen to verify (^0|) , so that the A = result is identical 
in the two cases. However, the effect of the correlation 
length is different for the two models: whereas for model 
A A retards both the NIOT and the NIDT, for model B 
the NIOT is retarded, but the NIDT is advanced. This 
is an indication of the nontrivial influence of the noise 
correlation length in the conserved case. 



V. STRONG COUPLING LIMIT 

In the limit of strong coupling, £) — > 00, the predic- 
tions of mean-field theory can be evaluated analytically 
and should agree with the results given by a standard 
linear stability analysis of the model. In order to ver- 
ify this agreement, we will now compute this limit for 
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the mean-field results obtained so far, for both models 
A and B. 



steepest-descent method, as before. Following the proce- 
dure described above for model A, we formally integrate 
now Eq. to obtain an equation for h, 



A. Model A 

In the mean-field approximation and in the limit D — > 
oo, the stationary probability distribution Pst{(f>, {4>)st) 
(pO|) becomes 



(0),t)=^((/'-(0)st), 



(41) 



as can be easily seen by means of a steepest-descent 
calculation. This expression verifies trivially the self- 
consistency relation (|2l]), which can thus no longer be 
used to determine (0)st- In order to do that, we now 
integrate Eq. ( p^ ) with respect to (j), and obtain 



{fmst + 'j'co (ff'((/.)<?(</.))st = 0. 



(42) 



For D — > oo, these averages are evaluated trivially using 
expression (^), and Eq. (^ ) becomes 



/(((/.),t,) + a2 co.g'((0).t)5((0)st) = O, 



(43) 



from which (0)st can be found. For model A and in the 
case defined by Eqs. the solutions of this equa- 

tion are either 



or 



1,2 



(<^)st = 



Co . 



(44) 



(45) 



This second set of solutions can only exist for a > —cr^ Cq. 
Hence, the critical point is given in this case by 



-o-^ Co : 



(46) 



in such a way that the ordered state appears for a > ac- 
The shift of the critical point increases the ordered re- 
gion, due to the effective multiplicative noise intensity 
tj^co. This shift, as seen in the previous Sections, in- 
creases with increasing noise intensity and decreases for 
increasing correlation lengths. This result coincides with 
the one given by a linear stability analysis of the homoge- 
neous state 16, 1^], as expected. However, in contrast 
with the linear stability analysis, this calculation can be 
extended to other situations and models not necessarily 
controlled by the linear term. 



B. Model B 

In this case, the stationary probability distribution 
( p6| ) given by the mean-field approach for each phase 
and for oo is also (E^) as can be seen using the 



h = -(/(0))st + 2rfa2(ci g(((/.),t) {g'mst - 

-co{g{(b)g'icj,)U). (47) 

At the limit _D — > oo, the averages appearing in the previ- 
ous expression are calculated using the stationary prob- 
ability distribution obtained above, leading to 

h = -/((0).t) + 2da^ (ci - Co) ff'((0)st) 5((0)st) . (48) 



In the case a < qt, the field is homogeneous and we 
can replace in the above expression {4))st — 4>o- Thus this 
equation gives us the value of h in this case as a function 
of the initial condition. The results for model B and in 
the case defined by (Q)-(0) are plotted in Fig. |^ versus 
the control parameter a, along with the values of h given 
by mean-field theory for finite but large D, obtained nu- 
merically in the previous Section. We can see that these 
mean-field results approach Eq. (^8|) as D increases, as it 
should be. The shift of the transition point increases for 
increasing coupling strength as can be seen from Fig. |^. 

We now turn to the case a > ut where h — 0. Now 
Eq. (|8|) can be solved for (0)st, which gives the values 
of the two bulk phases, 



(49) 



again for the particular model (0)-(0). The transition 
line is determined by setting — 0o hi the previous 
expression, which leads to 



ot = 00 ~ 2(icr^(co — ci) 



(50) 



and the critical point (for (f>Q = 0) is then 

Oc = — 2(itT^(co — ci) . (51) 

This result coincides with that coming from linear sta- 
bility analysis |16|. As in model A and in the previous 
Sections, the shift is in the direction of increasing the or- 
dered region. Due to the factor 2d this shift is larger than 
the one produced in model A for the same noise intensi- 
ties. Contrary to model A in the colored case, the shift 
does not depend only on the effective multiplicative noise 
intensity cr^ co but also on the noise correlation between 
first neighbors ci, which indicates the nontrivial influ- 
ence of the spatial correlation of the noise on conserved 
dynamics, as opposed to nonconserved dynamics where 
this influence disappears in the mean-field approach. 
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VI. CONCLUSIONS 



Mean-field theory has been previously applied to non- 
conserved models with additive and multiplicative white 
noises 22 . Here we have applied it in the case of 
spatially correlated multiplicative noise. Our mean-field 
results and numerical simulations of the complete model 
in two dimensions indicate the decrease of the ordering 
role of multiplicative noise when its correlation length 
increases. 

We have also extended mean-field theory to deal with 
conserved models by using the concept of a constant ef- 
fective field. As in the case of nonconserved systems, 
we have found that additive noise has a disordering role, 
whereas multiplicative noise has an ordering one. The 
latter increases for increasing multiplicative noise inten- 
sity and for decreasing noise correlation length. However, 
the quantitative effects of multiplicative noise are differ- 
ent in each model, the transition to order occurs earlier 
for model B than for model A. Moreover, mean field cal- 
culations show that the correlation length of multiplica- 
tive noise has nontrivial effects in the conserved case, 
while for model A it just decreases the effective noise 
intensity. Numerical simulations of the complete con- 
served model in two dimensions are in good agreement 
with mean-field predictions. 

Previous works on model A with additive and multi- 
plicative white noises have shown the presence of NIOTs 
and NIDTs. We have seen that, at least in the mean-field 
approach, these transitions appear for higher values of 
the noise intensity when multiplicative noise is spatially 
correlated. This is explained by the fact that the effective 
noise intensity decreases. Model B has also been found 
to go firstly through a NIOT and after through a NIDT 
when the multiplicative noise intensity is increased. As 
in model A, the NIOT is retarded when the correlation 
length of multiplicative noise increases. However, con- 
trary to what happens in model A, the NIDT is advanced 
what shows clearly different effects of noise correlation 
length upon conserved and nonconserved models. 

Finally, in the strong-coupling limit we have found an- 
alytical expressions for the critical-point shift and the 
steady-state bulk order parameter for both models A and 
B with additive white and multiplicative colored noises. 
These results coincide with previously reported predic- 
tions coming from linear stability analysis ||l6|| . 
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FIG. 1. Steady-state order parameter {4>)st versus control parameter a for model A. Lines are mean-field results and points 
correspond to numerical simulations of the complete model for a two-dimensional square lattice with 64 x 64 cells, mesh size 
Ax = 1, for system parameters X = 0.0 (circles and solid line), A = 0.5 (squares and dotted line) and X = 1.5 (triangles 
and dashed line). Other parameters are D = 3.7, e — 0.1 and = 5.0. All simulations in this paper use the same lattice 
parameters. 

FIG. 2. Steady-state bulk order parameter {(p) versus control parameter a for model B. Lines are mean-field results and 
points correspond to numerical simulations for X = 0.0 (circles and solid line), X — 0.5 (squares and dotted line) and X = 1.5 
(triangles and dashed line). Other parameters are D — 3.7, e — 0.1 and a'^ — 1.25. 

FIG. 3. Constant effective field h as a function of the control parameter a, as obtained from mean-field theory, for D = 3.7 
and (po ~ 0.2. All the cases with multiplicative noise have also additive noise with e — 0.1. 

FIG. 4. Steady-state bulk order parameters {</>) versus control parameter a for model A and B with additive white and 
multiplicative colored noises for different correlation lengths. Lines are mean-field results for model A with X = 0.0 (solid line) 
and X = 1.5 (dotted line) and for model B with X = 0.0 (solid line) and X — 1.5 (dashed line). Points correspond to numerical 
simulations of model A (empty symbols) and model B (full symbols). Of these, circles correspond to white multiplicative noise 
and triangles to X = 1.5. Other parameters are D = 3.7, e — 0.1, a\ = 5 and a% = 1.25. 

FIG. 5. Mean field steady-state bulk order parameters (0) versus multiplicative noise intensity for models A with X = 0.0 
(solid line) and X = 0.5 (dotted line) and for model B with X — 0.0 (solid line) and X = 0.5 (dashed line). Parameters are 
e = 1, a = 0.75 and D = 2.66. 

FIG. 6. Constant effective field h versus control parameter a for D cxi (solid line) as given by Eq. ^4^, (^i^d for D — 3.7 
(dashed line) and D = 20 (dotted line) coming from the mean-field approach described in Section IV. B. Other parameters are 
<t>o = 0.2, € = 0.1, (7^ = 1.25 and X = 0.5. 
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